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Abstract: A phenomenological extension of the well-known brane-world cosmology of Dvali, 
Gabadadze and Porrati (eDGP) has recently been proposed. In this model, a cosmological- 
constant-like term is explicitly present as a non- vanishing tension a on the brane, and an 
extra parameter a tunes the cross-over scale r c , the scale at which higher dimensional gravity 
effects become non negligible. Since the Hubble parameter in this cosmology reproduces the 
same ACDM expansion history, we study how upcoming weak lensing surveys, such as Euclid 
and DES (Dark Energy Survey), can confirm or rule out this class of models. We perform 
Markov Chain Monte Carlo simulations to determine the parameters of the model, using Type 
la Supernovae, H{z) data, Gamma Ray Bursts and Baryon Acoustic Oscillations. We also fit 
the power spectrum of the temperature anisotropies of the Cosmic Microwave Background 
to obtain the correct normalisation for the density perturbation power spectrum. Then, we 
compute the matter and the cosmic shear power spectra, both in the linear and non-linear 
regimes. The latter is calculated with the two different approaches of Hu and Sawicki (2007) 
(HS) and Khoury and Wyman (2009) (KW). With the eDGP parameters coming from the 
Markov Chains, KW reproduces the ACDM matter power spectrum at both linear and non- 
linear scales and the ACDM and eDGP shear signals are degenerate. This result does not 
hold with the HS prescription: Euclid can distinguish the eDGP model from ACDM because 
their expected power spectra roughly differ by the 3cr uncertainty in the angular scale range 
700 < I < 3000; on the contrary, the two models differ at most by the la uncertainty over 
the range 500 < £ < 3000 in the DES experiment and they are virtually indistinguishable. 
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1. Introduction 

The current ACDM concordance cosmological model excellently reproduces a large number 
of cosmological datasets. However, in order to explain the dynamics of galaxies and galaxy 
clusters and the Large-Scale Structures (LSS) of the Universe [1, 2, 3, 4, 5, 6], the observed 
pattern of the temperature anisotropies of the Cosmic Microwave Background (CMB) [7, 8] 
and the present accelerated expansion of the cosmos [9, 10, 11, 6, 12], we need to postulate 
that more than 80% of the matter in the Universe is made of weakly-interacting non-baryonic 
Dark Matter (DM) and that ~ 70% of the total energy budget is due to a particular negative- 
pressure Dark Energy (DE) in the form of a cosmological constant A. 

However, if one interprets the cosmological constant as the energy density of the vacuum 
of a quantum field, its measured value is more than hundred orders of magnitude smaller than 
what is predicted by quantum field theory (for a different perspective, see [13]). To avoid this 
problem, a number of alternative cosmological models has arisen over the last decades. They 
mainly try to reproduce the current data without adding any DM and/or DE, but rather by 
modifying or generalizing the law of gravity. 

Here, we focus on a particular class of these theories, known as brane-world cosmologies. 
They consider the Universe as a differential manifold, called bulk, with more than four di- 
mensions, which contains a 4D submanifold, the brane, with three spatial and one temporal 
dimensions; this submanifold is our spacetime. The model proposed by Dvali, Gabadadze 
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and Porrati [14, 15] is of particular interest, as its 5D Minkowskian bulk is able to induce 
a late-time acceleration of the Universe expansion on the brane in agreement with current 
observations. However, the growth of cosmological structures in the DGP model is strongly- 
suppressed compared to the ACDM model [16, 17, 18, 19, 20], and this also reflects on the 
DGP weak lensing power spectrum [21, 22]. 

A phenomenological extension of the DGP model has recently been proposed [23, 24, 25, 
26, 27]. In this extension, the DGP cross-over length r c , that determines the scale at which 
higher-dimensional gravitational effects become important, is tuned by a new parameter a, 
and a cosmological-constant term is explicitly present as a non-vanishing tension a on the 
brane. These models are free of ghost-like instabilities and one can obtain such a model with 
the so-called cascading gravity framework [28, 29, 30]. In this construction, our 4D brane lies 
within a succession of higher-dimensional DGP branes, embedded in each other within a flat 
bulk spacetime. 

The eDGP model we use has been proposed as a generalisation of the DGP model in 
a spatially flat Universe, but it might be also worth investigating the non-flat case. For 
the standard DGP, which presents strong tension between SNela data and CMB distance 
indicators, open models provide a better fit to the geometric data, because the distance to 
the last scattering surface is not constrained by the flat geometry. By introducing VLk as an 
additional free parameter, the constraints on the cosmological parameters would be weaker 
and, consequently, the importance of the tests on the growth of structures would increase. 
Unfortunately, DGP generalisations for the non-flat case similar or related to the model we 
examine here still lack. 

In the present work, we explore this extended DGP (eDGP) model in detail. We perform 
Monte Carlo Markov Chains to obtain the best-fit value of the cosmological parameters. To 
this end, we use a number of observational datasets, which have been shown to be useful in the 
case of modified gravity models [31]. Specifically, we use: the Hubble constant Hq, estimated 
by the HST Project using a well calibrated set of local distances [32]; the latest 397 Type 
la Supernova (SNIa) distance moduli of the Constitution sample, compiled in Ref. [33]; the 
updated Gamma Ray Burst (GRB) Hubble diagram recently presented in Ref. [34]; Baryon 
Acoustic Oscillation (BAO) data from the seventh release of the Sloan Digital Sky Survey 
(SDSS); and distance priors from WMAP5. We also run the chains for the ACDM model, to 
compare our results. We also perform a chi-squared analysis on the temperature anisotropy 
power spectrum of the CMB with the most recent data coming from WMAP7 [12], to estimate 
the normalisation of the power spectrum of the matter density perturbations. 

We calculate the matter power spectrum in both the linear and the non-linear regime of 
perturbations, the latter with both the recent fitting formulae of Khoury and Wyman [35], 
who obtained a useful modification of the HALOFIT procedure [36] for the eDGP model, and 
the interpolation proposed by Hu and Sawicki [37]. 

Finally, we compute the weak lensing convergence power spectrum. Weak lensing is a 
powerful tool [38, 39, 40, 41, 42], and it is particularly useful to aim of testing and constraining 
modified/alternative cosmological models [43, 44, 22, 45, 46, 47, 48]. Gravitational lens effects 
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are due to the deflection of light occurring when photons travel near matter, i.e. in the 
presence of a non-negligible gravitational field. Therefore, the cosmic convergence and shear 
encapsulate information about both the sources that emit the light and the structures that 
the photons cross before arriving at the telescope: weak lensing enables us to explore both 
the basis of the cosmological model and the LSS of the Universe. In other words, it provides 
information about geometry and dynamics. The study of the power spectrum of weak lensing 
can be a crucial test, particularly with a view to upcoming weak lensing surveys such as 
the Dark Energy Survey (DES) 1 [49, 50], the Panoramic Survey Telescope Rapid Response 
System (Pan-STARRS) 2 [51, 52] and future space surveys such as Euclid 3 [53, 54]. 

The paper is organized as follows: in § 2 the eDGP model is outlined; in § 3 we present the 
cosmological parameters obtained with Monte Carlo Markov Chains and the normalisation 
of the power spectrum of the density perturbations derived from the fit to the temperature 
anisotropies of the CMB; in § 4 we show our results: the matter power spectrum (§ 4.1) and 
the weak lensing shear (§ 4.2). Finally, in § 5 conclusions are drawn. 

2. The extended DGP model 

In the DGP model [14, 15, 55, 18], the bulk has five dimensions and is Minkowskian, and 
the brane is flat. In this case, the more general homogeneous and isotropic metric has three 
different scale factors, one for the time t, one for the three spatial coordinates x = (r, <p) 
and one for the fifth dimension y, and takes the form 4 

9AB dx A dx B = -N 2 (r, y) dt 2 + A 2 (r, y) dr 2 + B 2 (r, y) dU 2 + dy 2 , (2.1) 

where dO 2 = r 2 dt? 2 + r 2 sin 2 $ dip 2 is the differential solid angle. 

In this model, gravity is general relativity (GR), but the Hilbert-Einstein Lagrangian is 
given in five dimensions. Thus the coupling constant of the gravitational interaction is no 
longer the Planck mass Mp, because gravity is spread into an additional dimension. This new 

M 2 

energy scale M^ + i sets a cross-over length r c = 2M ^a that represents the scale at which the 
true, five dimensional gravitational effects become dominant with respect to the Einsteinian 
four dimensional gravity. 

The model is constrained by putting all the matter and energy species on the brane, and 
by letting the fifth dimension be filled only by gravity; the Einstein equations thus give a 
different relation for the evolution of the scale factor a(t) = A(t, y = 0) on the brane. The 
standard Hubble parameter is H = a/a, where a dot denotes a derivative with respect to the 

x http: / /www. darkenergysurvey.org 
2 http://pan-starrs. ifa.hawaii.edu 

3 http://sci.esa.int/science-e/www/area/index.cfm?fareaid=102 

4 We use units such that c = 1 and uppercase Latin indices run over bulk dimensions, Greek indices over 
4D spacetime dimensions, whereas lowercase Latin indices label spatial coordinates. 
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proper time, and the expansion history as a function of the redshift z reads 

2 



+ Q K (l + z)\ (2.2) 



where i = K,m,r stands for curvature, matter or radiation, with rij = 2,3,4, respectively, 
and Q a = 8irGa/ (3i?o 2 ) is a cosmological-constant term which can be interpreted as a non- 
vanishing tension a of the brane. We expect that, in the higher-dimensional framework, Q, a 
could find a natural origin. 

Observations indicate that the Universe is spatially flat [5]; therefore, the DGP model 
has been deeply investigated in the case of Qk = 1 — $7 m — Q r — Q, a — l/(r c Ho) = 0. In this 
case, the modified Friedmann equation is 

H 2 ±y = H 2 ^«~ n< + • (2.3) 

Hereafter, we will refer to the upper sign in the modified Friedmann equation as the "normal" 
branch, and to the lower sign as the "self-accelerating" branch. Indeed, with the choice of 
the minus sign, the DGP model can achieve an accelerated expansion without a cosmological 
constant, unlike the case of the normal branch, where a non-null f2 CT term is needed. 

In this paper we focus on a phenomenological extension of DGP model [23, 24, 25, 26, 27]. 
Here, the modified Friedmann equation (2.3) is generalised as 

H 2 ± = Ho 2 + Ox J , (2.4) 

where a is a free parameter which is strictly related to the graviton propagator [23, 26]. The 
condition Q.k = clearly implies 

n a = i - n m - n r ± (r CJ ff ) 2(a_1) . (2.5) 



Such a model presents a graviton propagator which is proportional to 

1 



k 2 + k 2a r c 2 ^ a - 1 ) 



(2.6) 



This would follow from a phenomenological model of modified gravity proposed in [24] and 
it is a power-law generalisation of the graviton propagator in the DGP brane-world. Since 
we are interested in a long distance modification of gravity, we assume that a < 1. The 
unitarity constraint requires a > [24]. It can be seen that a = 1 can be absorbed into a 
redefinition of the Newtonian gravitational constant G. The propagator with a -C 1 has some 
connection to the so-called cascading DGP brane-world [28] and thus the Friedmann equation 
with a <C 1 might be obtained in such higher codimension models. However, we would like to 
stress that a detailed analysis of higher codimension DGP models has yet to be undertaken 
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and no brane- world model has been known so far that leads to Eq. (2.4). Therefore, we shall 
view Eq. (2.4) as a phenomenological starting point for our modified gravity. It is easy to 
see that the standard DGP model is recovered when a = 1/2, while the choice a = 0, often 
called "degravitation," leads to an expansion history identical to that of the ACDM model. 



3. Parameter estimation 



We now estimate the model parameters required to reproduce the observed expansion history 
of the Universe. In the Bayesian approach to model testing, we explore the parameter space 
through the posterior probability density function 

7T(/l)L S NIa(p)^GRB(p)^BAo(p)-^WMAP5(p), (3.1) 

where p denotes the set of model parameters, ir(h) is the prior on the Hubble constant h and 
Ld is the likelihood function related to the dataset D. 
The Gaussian prior 

i ( h nsT- h \ 2 

ir(h) oc e 2 V -hst J (3.2) 

stems from the results of the HST Project [32] which estimated the Hubble constant Hq using 
a well calibrated set of local distances. Averaging over the different methods, the survey gives: 

/ihst ± o-hst = 0.72 ± 0.08 (3.3) 

which is independent of the cosmological model. 

The second and the third terms in the likelihood (3.1) are both related to the Hubble 
diagram, and involve SNela and GRBs. The SNela provided the first piece of evidence of the 
cosmic acceleration based on the predicted distance modulus 

Mth(*|p) = 25 + 5 log [(1 + z) X (z\p)} , (3.4) 

where dx = dz/H(z) is the differential radial comoving distance. The likelihood function is 
then defined as 

SNIa ^ \J det (27rC S Nia) ' ^ ^ 

where A/x is sl A/sNia~ vector containing the values /iobs 

(z^ - Hth{zi) and CsNia is the A/sNia x 
A/sNia covariance matrix of the SNela data. Here, A/sNia is the total number of SNela used. 

SNela are limited to z ~ 1.5. Therefore we have to resort to a different distance indicator 
to probe the Hubble diagram to higher redshift. Thanks to the enormous energy release that 
makes them visibile up to z ~ 6.6, GRBs stand out as possible candidates for this purpose. 
As a result, Schafer [56] has provided the first GRB Hubble diagram containing 69 objects 
with n hs(z) estimated by averaging over 5 different 2D correlations. Here we use the updated 
GRB Hubble diagram recently presented in [34] based on a model-independent recalibration 
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of the same 2D correlations used by Schafer. Since there is no correlation between the errors 
of different GRBs, the likelihood function simply reads 



_ 1 V %RB 
2 Z-ii=l 

^grb(p) oc e 



Mobs( z i)-Mth( z tlP) 
V /<T i +CT GRB 



(3.6) 



where A/grb is the total number of GRBs and cgrb quantifies the intrinsic scatter inherited 
from the scatter of GRBs around the 2D correlations used to derive the individual distance 
moduli. 

Lbao(p) is the likelihood function related to the BAO acoustic scale, which is effectively 
constrained by the detected peak in the correlation function of the luminous red galaxies. If 
we define the effective distance [57, 58] 



d v (z) 



X 2 (z)- 



(3.7) 



the ratio d z = r s (zd)/dv(z), where r s (z) is the sound horizon at the drag epoch Zd, is a 
well constrained quantity for z = 0.2 and z = 0.35 and its values are do.2 = 0.1905 and 
d .3 5 = 0.1097 [57, 58, 8]. 

Finally, the last term in the likelihood (3.1) is related to the properties of the CMB. 
Following [8], we define three fitting parameters for comparison with the WMAP5 data, i.e. 
the redshift of the surface of last scattering z rec , the shift parameter at that epoch 

-R(^rcc) = V / ^m- H o(l + 2 re c)<U(-Zrec), (3.8) 

and the acoustic scale 

l a = (l + z rcc )^j^. (3.9) 

Here, d^iz) = x( z )/(^ + z) is the proper angular diameter distance. We compare our results 
with the values z rec = 1090.04, R(z rec ) = 1.71 and l a = 302.1. 

The constraints on the individual parameters of the ACDM, eDGP normal and self- 
accelerating branches (eDGPn and eDGPs, respectively) are summarized in Tables 1-3, where 
we give the best-fit, mean and median values and the 68% and 95% confidence limits. These 
constraints are obtained for the entire set of cosmological observational data we choose for 
this work, i.e. SNela, H(z) data, BAOs and GRBs. 

We find that the eDGP model is able to fit these datasets. Moreover, in Fig. 1 we show 
how the expansion rates in ACDM and in the eDGP model differ by less than 0.001% over a 
wide range of redshift. 

3.1 CMB anisotropy power spectrum 

The CMB probes the geometry of the background expansion as well as the formation of LSS. 
To predict these effects we use the parameterised-post Friedmannian [37] modifications to the 
CAMB (Code for Anisotropies in the Microwave Background) routines 5 [59, 60, 61]. 

5 http : //camb . info/ 
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ACDM 



X 


fit 


(x) 


•^med. 


68% CL 


95% CL 




0.28 


0.28 


0.28 


(0.27, 0.29) 


(0.255, 0.31) 


io 4 n r 


1.1 


1.1 


1.1 


(1.0, 1.1) 


(9.9, 1.1) 


h 


0.71 


0.71 


0.71 


(0.69, 0.72) 


(0.67, 0.74) 


cgrb 


0.40 


0.37 


0.38 


(0.27, 0.47) 


(0.10, 0.57) 



Table 1: Constraints from SNela, H(z), GRBs and BAOs on the ACDM model parameters. 

eDGPn 



X 


S'bcst-fit 


(x) 


3?mcd. 


68% CL 


95% CL 




0.28 


0.28 


0.28 


(0.27, 0.29) 


(0.25, 0.31) 


10 4 f) r 


1.1 


1.1 


1.1 


(1.0, 1.1) 


(9.9, 1.1) 


r c H 


52.375 


155.04 


87.16 


(28.20, 284.94) 


(5.72, 766.268) 


a 


0.098 


0.116 


0.083 


(0.033, 0.187) 


(0.013, 0.481) 


h 


0.71 


0.71 


0.71 


(0.69, 0.72) 


(0.67, 0.74) 


cgrb 


0.42 


0.375 


0.38 


(0.28, 0.48) 


(0.18, 0.56) 



Table 2: Constraints from SNela, H(z), GRBs and BAOs on the eDGPn model parameters. 

eDGPs 



X 


■^best-fit 


(x) 


•^med. 


68% CL 


95% CL 




0.28 


0.28 


0.28 


(0.27, 0.295) 


(0.25, 0.31) 


io 4 n r 


1.1 


1.1 


1.1 


(1.0, 1.1) 


(9.9, 1.15) 


r c H 


19.43 


88.245 


73.33 


(24.20, 146.10) 


(3.225, 287.58) 


a 


0.517 


0.103 


0.077 


(0.029, 0.166) 


(0.013, 0.400) 


h 


0.71 


0.71 


0.71 


(0.69, 0.72) 


(0.67, 0.74) 


CGRB 


0.415 


0.36 


0.37 


(0.25, 0.48) 


(0.90, 0.56) 



Table 3: Constraints from SNela, H(z), GRBs and BAOs on the eDGPs model parameters. 



We perform a chi-squared analysis with the most recent data coming from WMAP7 [12] 
to get the best-fit value for the temperature anisotropy power spectrum normalisation, for 
the ACDM model and both branches of the eDGP model. We thus obtain the normalisation 
of the matter power spectrum that we use in the next section and in further calculations. The 
three fits of our models to the CMB data are indistinguishable: with the WMAP7 best-fit 
value n s = 0.963 for the tilt of the primordial power spectrum, the three reduced chi-squared 
are XeDGPn = XeDGPs = Xacdm = 1-157. We derive ag = 0.858, the rms mass fluctuations on 
the scale R = 8/i _1 Mpc. 
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z 

Figure 1: Relative error between the ACDM and the eDGP expansion histories. 

4. Power spectra 

4.1 Matter power spectrum 

A faster growth of structure at early times, followed by a slower structure growth from 
the point on where the modification of gravity becomes important, is a feature common to 
eDGP models and several modified gravity theories. As usual, we define the growth factor 
D + (a) = 5(x, a)/<5(x, a = 1), where 5 = 5p/p is the density contrast. The evolution equation 
of D + in the matter-dominated epoch for DGP brane-worlds becomes [16, 62, 18, 19, 63, 25] 

6 + 2H5 = ^n m H 2 (l + A, (4.1) 

where 5 = Sp/p is the density contrast and 

13 = 1 ± 2 (r c H) 2(1 - a) 



10 

10 
10 



11 

12 



i + 



2(1 - a)H 
3IP 



(4.2) 



We have shown that the Hubble parameter in both branches of the eDGP model differs 
by less than 0.001% from that of the ACDM model (Fig. 1), and the normalisation of the 



primordial power spectrum obtained from the CMB temperature anisotropy analysis is again 
the same in both eDGPs and eDGPn (§ 3.1). Therefore, hereafter our results will refer only 
to the eDGP model, because the curves in the two branches are almost indistinguishable. 

4.1.1 Linear regime 

By performing the Fourier transform <5^(a) of the density contrast, we can construct the 
present-day matter power spectrum P(k) = P s (k,a = 1), where, for a generic field fk(a) 



with 5d the Dirac 5 function. In the DGP model, the Newtonian potential 3>fc(a) an d the 
metric potential *fc(a) depend on the scale k and thus their evolution deeply differs from 
the corresponding quantities in the ACDM model; however, the growth factor does evolve 
similarly to ACDM for a wide range of scales [20]. Thanks to this peculiarity, we can use the 
growth factor D + (a), which is solution of Eq. (4.1), to construct the matter power spectrum 
as 



Here, n s is the tilt of the primordial power spectrum, 5h is its normalisation and T{k) is the 
matter transfer function which describes the evolution of perturbations through the epochs 
of horizon crossing and radiation-matter transition. As mentioned in § 3.1, we use n s = 0.963 
and 8h such that as = 0.858. The transfer function is in Ref. [64]. 

4.1.2 Non-linear regime 

To obtain the non-linear matter power spectrum, whose contribution is particularly relevant 
for computing the weak lensing cosmic shear signal, we use two different prescriptions. On 
one hand, A r -body simulations in eDGP models and degravitation theories have recently 
been performed [35] and they show that, in order to correctly reproduce the matter power 
spectrum, the halofit fitting formula? [36] have to be slightly modified: the new set of 
parameters works well for the whole range of a, although additional A-body simulations 
covering a wider volume in the parameter space still lack. The new parameters of the fitting 
formulae are listed in Table 4. The parameter 3 is related to a redefinition of the y parameter 
of halofit, i.e. y — > <5y. We will denote this non-linear perscription as "KW" from Khoury 
and Wyman (2009), Ref. [35]. 

For modified gravity to agree with Solar system observations, the non-linear matter power 
spectrum has to approach the standard ACDM solution on small scales. This means that 
the non-linear power spectrum has to be an interpolation of two power spectra. The first 
is the modified gravity non-linear power spectrum Py lG (k,z), which is obtained without the 
non-linear interactions that are responsible for the recovery of GR. This is equivalent to 
assume that gravity is modified down to small scales in the same way as in the linear regime. 
The second term, Pq R (A;, z), is the non-linear power spectrum obtained in the DE model 



(fk(a)f k >*(a)) = (2vr) 3 5 D (k - k') pf(k, a), 



(4.3) 




(4.4) 



9 



log 10 a n 0.841og 10 < td 
log 10 6 n log 10 6^ td + log 10 l.l 
log 10 c n log 10 4 td + log 10 1.05 

logio logic /"n d + lo gl0 °- 875 

log 10 u n log 10 + log 10 0.875 
a n 0.8< td 

S 1.035 

Table 4: Parameters for the modified non-linear fitting algorithm from [35] (left column) as a function 
of the parameters in the standard HALOFIT formula [36] (right column) . 



that follows the same expansion history of the Universe as the modified gravity model, yet 
obeying to GR. In other words, this is the non-linear P s (k, z) which will have a ACDM model 
with an expansions history H(z) equivalent to that of the modified gravity theory. Since the 
Hubble parameter in the eDGP model with the parameter values obtained with our Markov 
Chains is the same as that in the ACDM model, Pq R (&, z) is simply the ACDM matter power 
spectrum. We use the fitting formula proposed by Hu and Sawicki (2007), Ref. [37] 

pS (h Pj iG {k,z) + c nl {z)T?{k,z)P^{k,z) 

where T?{k,z) = [k 2 P s {k, z)/2vr 2 ] ai picks out non-linear scales, since P s (k , z) is the linear 
power spectrum of Eq. (4.3); c n \(z) = ^4(1 + z) 0-2 determines the scale at which the power 
spectrum approaches the GR result as a function of redshift. Their functional forms have 
been obtained by perturbation theory [65] and confirmed by A^-body simulations [66, 67]. 
Here, we use a\ = 1, 02 = 0.16 and A = 0.3. We will denote this nondinear perscription as 
"HS." 

In addition to the two non-linear prescriptions mentioned above, Schmidt et al. [68] 
have recently studied the spherical collapse model in the DGP brane-world gravity. They 
show how results on the mass function, halo bias and power spectrum obtained with N- 
body simulations [67, 69] can be described semi-analytically. However, they find a tension 
between the simulations and their reconstructed power spectrum. Specifically, the P n \(k) in 
the self-accelerating branch matches the simulated data very well, whereas there is a 10-^30% 
discrepancy in the normal branch. Nonetheless, their spherical collapse predictions appear to 
be able to provide an effective way of including non-linear effects in results obtained from DGP 
simulations in the linear regime. Further investigations are worth performing, in particular 
in the perspective of generalising such an approach to the eDGP model. 

In Fig. 2 we show the non-linear matter power spectra for the ACDM and the eDGP 
models. For the latter, we show the power spectra obtained with both the KW and the HS 
methods described above. As expected, the eDGP-HS non-linear P(k) approaches the ACDM 
signal at large A;'s, because the prescription was explicitly proposed to track the non-linear 
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behaviour of GR at small scales. Unexpectedly, we find that the KW prescription also is 
able to reproduce the non-linear ACDM power spectrum. This result was not found in [35], 
because our a and r c , unlike theirs, derive from a Markov Chain Monte Carlo analysis which 
is based on the data describing the expansion history of the real Universe. 
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Figure 2: Ratio of the non-linear matter power spectra of the ACDM and the eDGP models. For 
the eDGP model, the non-linear power spectra obtained with [37, 35] are shown. 

4.2 Weak lensing signal 

It is known from GR that light beam paths are curved by the presence of matter. In the 
weak lensing framework the deflection of light is small and, consequently, we can use Born's 
approximation, where lensing effects are evaluated on the null-geodesic of the unperturbed 
(unlensed) photon [38, 70]. 

In the DGP model, and similarly in its phenomenological extensions, for the two poten- 
tials |$| |^| holds. Thus, the relation between the distribution of matter overdensities in 
the Universe and the two metric perturbations, the potentials <3? and is not trivial. In 
GR, in the matter-dominated era, when there is no anisotropic stress, <3? = — and therefore 
we can simply use the Newtonian potential <I>, thanks to the canonical Poisson equation, to 
compute cosmic convergence and shear. However, in general, the weak lensing effect is due 
to the combination of both the Newtonian and the metric potential. We will refer to this 



- II - 



combination as the "deflecting potential," and we will denote it with 6 

$ _ \Jr 

T ee — . (4.6) 

A technique for obtaining the Fourier modes of the deflecting potential starting from the 

growth factor in the extended DGP model has recently been proposed [35]. This modified 
Poisson equation reads 

T fc (o) - -n m H — p— - , (4.7) 

where 

»(*,.)_ G 2( izlMBzi) ,4.8, 

and 

e(fe,o) =8 5 2 #Va n ,5 fc , (4.9) 

. = ~ (4.10) 

P 5 (/c, a) is the matter power spectrum as a function of the comoving wavenumber k and the 
scale factor a, as computed in § 4.1. 

All weak lensing observables may be expressed in terms of the projected potential 

<Kn) = |d X ^T(n, X ) (4.11) 

where 

W(X) = -X d X ' ^^n(x') (4.12) 

is the weight function of weak lensing, with n [xi z )} representing the redshift distribution of 
the sources, such that Jd%n(x) = 1- We now introduce a distortion tensor [71, 70] 

<^(x) = /VxWJT^n.x'), (4.13) 

JO 

where commas denote derivatives with respect to directions perpendicular to the line of sight. 
The trace of the distortion tensor represents the cosmic convergence 

/e(x) = £ii(x) + 0,22(x) (4.14) 

and the (complex) shear is the linear combination 

7 (x) = 0,n (x) - 0, 22 (x) + 2#,i 2 (x). (4.15) 



6 In the literature, the combinations of the two potentials are often indicated with $± = — ($ ± *&) /2. 
However, different authors use to refer to the metric potentials in the length element differently. For the sake 
of simplicity, we have chosen this notation, because we are only interested in the combined effect responsible 
for the weak lensing signal. 
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In the flat-sky approach we expand the projected potential (f>(n) in its Fourier modes 



m = Jj^f met£h - (416) 
The power spectrum is defined as the Fourier transform of the 2D correlation function 

(</>*{£)(/>(£')) = (2it) 2 5 d (£ - £')C^(£). (4.17) 

Thus, we have [72] 

C**{t) = f°d X ^Mp T (-,x) (4.18) 

where P r (k, z) is the power spectrum of the deflecting potential constructed following Eq. (4.3), 
and we have introduced Limber's approximation, where the only Fourier modes that con- 
tribute to the integral are those with I » k\- The cosmic convergence and shear can be 
obtained with 

C KK (£) = C^(£) = tC H (t). (4.19) 

We compute the shear power spectrum of light emitted by background galaxies, because 
this effect is relevant for the study of alternative/modified gravity theories. It has been shown 
that, although the CMB convergence can give a signal similar to that of the ACDM model 
[45], the background galaxy power spectra can be different. This is due to the simple fact 
that models alternative to the standard ACDM model have to reproduce the same behaviour 
in the early Universe, where CMB constraints are stricter, but can differ at late times. 

We calculate results for realistic oncoming surveys: a ground-based survey similar to 
DES, and a space-based survey such as Euclid, using the redshift distributions shown in 
Fig. 3. The redshift distribution for our ground-based survey is chosen to be the same as for 
CFHTLS [73] giving a median redshift z m = 0.826, a sky coverage ©deg = 5,000deg 2 and a 
mean galaxy number density n = 13.3 arcmin~ 2 ; for Euclid we use the distribution given by 
[74, 75] with z rn = 0.91, @ 2 dcg = 20,000deg 2 and n = 35arcmin- 2 . 

Figs. 4-5 show the shear power spectra 1(1 + l)C 77 (£)/(27r) of light from background 
galaxies in the eDGP model and in the standard ACDM model. The shear signal shows the 
same behaviour of the non-linear matter power spectrum shown in Fig. 2: the weak lensing 
signal in eDGP model appears to be suppressed with respect to ACDM, when using the 
HS procedure for computing the non-linear P s (k,z), whereas the eDGP-KW shear signal is 
virtually indistinguishable from ACDM. 

We also show the expected errors as computed with [72, 71] 



where f s ^ y = 2 lcg 7r/129, 600 is the fraction of the sky covered by a survey of area and 
(Tint 2 ) ' 5 — 0.4 is the galaxy-intrinsic rms shear in one component. The errorbars clearly show 
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Figure 3: The source redshift distributions used: for ground-based surveys (DES) with median 
redshift z m = 0.826 and for a space-survey (Euclid) with median redshift z m = 0.91. 

that the depth and the wide field of a space-based survey is necessary to discriminate between 
the ACDM and the extension of the DGP model. In fact, DES could do so, in principle, only 
at Id level in the range of angular scales 500 < i < 3000, whereas, thanks to its wide and 
deep field, Euclid can achieve 3a level in the range 700 < £ < 3000. 

5. Conclusions 

We investigate a phenomenological extension of the Dvali-Gabadadze-Porrati (DGP) brane- 
world cosmology, which presents an extra parameter a as an exponent of the non-quadratic 
term in the modified Friedmann equation. The DGP model admits two Hubble parameter 
H(z) solutions of the modified Friedmann equations: one of them presents a self-accelerating 
behaviour at late-times even without any cosmological constant term, hence it is known as the 
"self-accelerating" branch. The other solution is often called the "normal" branch. Similarly, 
the extended DGP model has two expansion histories of the Universe, and we refer to them as 
eDGPn and eDGPs, for the normal and self-accelerating branch, respectively. In this work, 
we allow the eDGP model to have an explicit cosmological-constant-like term in the modified 
Friedmann equation, which can be interpreted as a non-vanishing tension a on the brane. 
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Figure 4: Shear power spectra £(£ + 1)C 77 (£) / (2tt) in the eDGP model and in the standard ACDM 
model of photons coming from background galaxies expected from a ground-based DES-like survey dis- 
tribution of sources with median redshift z m — 0.826: blue/dark with the KW non-linear prescription 
[35], green/light with HS [37]. Error bars are la. 

We compute Monte Carlo Markov Chains to estimate cosmological parameters consistent 
with current observations. We select a sample of cosmological observables. In particular, we 
use H(z) data, Type la Supernova? (SNela) and Gamma Ray Bursts (GRBs) as standard 
candles, Baryon Acoustic Oscillations (BAOs) from SDSS and distance priors from WMAP5 
data. We run the Markov Chains for both branches of the eDGP model and for the ACDM 
model, in order to make direct comparisons between the two models. We find that the eDGP 
models can fit current data at least as well as the standard ACDM model. 

By using an extension of the CAMB code with the parameterised-post Friedmannian 
framework for modified gravity, we also obtain the CMB temperature anisotropy power spec- 
trum for the eDGP models. With the cosmological parameters obtained with our Markov 
Chains, we find a perfect agreement between the recent WMAP7 data and the power spec- 
tra, for both ACDM and the eDGP models. Specifically, the chi-squared analysis gives 
XeDGPn = AeDGPs = Xacdm = 1-157. The only parameter left free by the Markov Chain 
Monte Carlo simulations is as, the rms mass fluctuations in a sphere with radius 8/i _1 Mpc. 
Our chi-squared analysis of the CMB power spectrum yields as = 0.858. 
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Figure 5: Shear power spectra £(£ + 1)C 77 (£) / (2tt) in the eDGP model and in the standard ACDM 
model of photons coming from background galaxies expected from a space-based Euclid-like survey 
distribution of sources with median redshift z m = 0.91: blue/dark with the KW non-linear prescription 
[35], green/light with HS [37]. Error bars are la. 

With this normalisation, we compute the matter power spectrum of the density pertur- 
bations in the linear regime. When we correct the linear power spectrum by using the recent 
modification to the halofit fitting formulas of Ref. [35] (KW) for the eDGP model, we see 
that, at non-linear scales, unlike the results reported in [35], the matter power spectrum in 
both branches of the eDGP model approaches the ACDM curve. This result is a consequence 
of the fact that our a and r c parameters derive from the Monte Carlo Markov Chains applied 
to the current data of the Universe expansion history. We also use the non-linear interpolation 
of Ref. [37] (HS), which has been explicitly proposed to track the non-linear behaviour of GR 
at small scales. 

Finally, we compute the weak lensing shear power spectrum as it is expected from a 
space-based mission such as Euclid and from a ground-based survey like DES, and we show 
that the eDGP model is difficult to distinguish from ACDM: indeed, at both large and very 
small angular scales the observational constraints are loose and the differences between the 
two models are within the uncertainties. Specifically, if the correct non-linear prescription 
were KW, the ACDM and eDGP shear signal would be degenerate. On the contrary, with 
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the HS non-linear prescription, the ACDM and the eDGP models differ. However, DES is 
unable to discriminate between the two models because they differ by at most la in the range 
of angular scales 500 < i < 3000. On the other hand, the depth and the wide field of a 
space-based survey such as Euclid can in principle discriminate between the eDGP and the 
ACDM models, because they differ by 3a in the range 700 < I < 3000. 
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